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LARGE AMOUNTS OF OPTICALLY-OBSCURED STAR FORMATION IN THE HOST GALAXIES OF SOME 

TYPE-2 QUASARS. 
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ABSTRACT 

We present Hubble Space Telescope images, and spectral energy distributions from optical to infrared 
wavelengths for a sample of six 0.3 < z < 0.8 type-2 quasars selected in the mid-infrared using data 
from the Spitzer Space Telescope. All the host galaxies show some signs of disturbance. Most seem 
to possess dusty, star-forming disks. The disk inclination, estimated from the axial ratio of the hosts, 
correlates with the depth of the silicate feature in the mid-infrared spectra, implying that at least 
some of the reddening towards the AGN arises in the host galaxy. The star formation rates in these 
objects, as inferred from the strengths of the PAH features and far-infrared continuum, range from 3 — 
9OM yr _1 , but are mostly much larger than those inferred from the [On] 3727 emission line luminosity, 
due to obscuration. Taken together with studies of type-2 quasar hosts from samples selected in the 
optical and X-ray, this is consistent with previous suggestions that two types of extinction processes 
operate within the type-2 quasar population, namely a component due to the dusty torus in the 
immediate environment of the AGN, and a more extended component due to a dusty, star forming 
disk. 

Subject headings: quasarsigeneral - galaxies :Seyfert - infrared:galaxies - galaxies: starburst 



1. INTRODUCTION 

Co-evolution of black holes and their host galaxies is 
a prediction of most models for the observed close cor- 
relation between black hole mass and bulge luminos- 
ity/velocity dispersion in local galaxies (e.g. Kauffmann 
& Haehnelt 2000; di Matteo, Springel & Herquist 2005). 
Studies of local AGN and type-1 (unobscured) quasars 
show a correlation between AGN luminosity in the mid- 
infrared and starburst luminosity, estimated from the 
far-infrared emission and/or PAH features (Schweitzer 
et al. 2006; Shi et al. 2007; Maiolino et al. 2007), and 
the host galaxies of quasars with strong far-infrared ex- 
cesses do seem to be forming stars (Canalizo & Stockton 
2001). However, in general, host galaxy imaging studies 
in the optical and near-IR do not suggest that strong star 
formation occurs simultaneously with quasar activity in 
the majority of quasars. Imaging of type-1 quasar host 
galaxies show that the hosts of luminous quasars tend 
to be elliptical galaxies with little obvious signs of dis- 
turbance in Hubble Space Telescope (HST) images (e.g. 
Floyd et al. 2004), and only about 30% of quasars show 
signs of mergers or interactions in HST or adaptive optics 
images (Marble et al. 2003; Guyon et al. 2007), though 
deeper (several orbit) HST/ ACS imaging often reveals 
faint tidal features in quasar hosts (Canalizo et al. 2007; 
Bennert et al., in prepraration). Studies of type-2 (ob- 
scured) quasars have so far been consistent with this pic- 
ture, and with the idea that the obscuration takes place 
in a nuclear torus. Sturm et al. (2006), in their Spitzer 
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study of X-ray selected type-2 quasars, make the claim 
that "type-2 quasars are not ULIRGs" based on feature- 
less mid-infrared spectra. HST imaging studies of the 
host galaxies of type-2 quasars selected from the Sloan 
Digital Sky Survey (SDSS) show that they are consis- 
tent with being hosted by relatively undisturbed ellipti- 
cal galaxies, once the likely effects of scattered light are 
accounted for (Zakamska et al. 2006). 

In this Letter, we present results of Spitzer spec- 
troscopy and photometry, and Hubble Space Telescope 
(HST) imaging of a sample of type-2 quasars, and dis- 
cuss the implications for the co-evolution of galaxies and 
black holes. Our sample is selected in the mid-infrared, 
which is less sensitive to reddening in the host galaxy 
than X-ray or optical selection, and their type-2 nature 
allows us to image the central regions of the host with 
high fidelity. The combination of infrared SEDs and op- 
tical imaging of this sample thus allows us for the first 
time to present a coherent picture of the nature of the 
star formation in the host galaxies of type-2 quasars. 

2. OBSERVATIONS 

We selected all six candidate type-2 quasars with 0.3 < 
z < 0.8 from the lmJy 8/im flux density limited sample 
of Lacy et al. (2004). The redshift range ensures that 
these objects are luminous enough to be on or above the 
quasar /Seyfert divide in the mid-infrared (5 jam luminos- 
ity, L 5fim ~ 10 23 - 6 WHz-\ Lacy et al. 2005a), but low 
enough for rest-frame £?-band to be within the F814W 
filter of the Advanced Camera for Surveys (ACS) aboard 
HST. The objects were observed with the Multiband 
Imaging Photometer (MIPS) and Infrared Spectrograph 
(IRS) as part of Spitzer program 20083, and with HST as 
program GO- 10848 (Table 1). Redshifts were obtained, 
and optical spectral classifications attempted (Lacy et al. 
2007). Low dispersion spectra were also obtained with 
the Spex instrument (Rayner et al. 2003) in prism mode 
at the Infrared Telescope Facility (IRTF). 



2 



Lacy et al. 



TABLE 1 

Observations 



Name 


IRS 


MIPS 


ACS filters 




reqkey 


reqkey 




SSTXFLS J172123. 1+601214 


14016768 


14018304 


F625W,F435W 


SSTXFLS J171147.4+585839 


14016000 


14017536 


F814W,F555W 


SSTXFLS J171324. 2+585549 


14016256 


14017792 


F775W,F475W 


SSTXFLS J171831. 5+595317 


14016512 


14018048 


F775W,F475W 


SSTXFLS J171106. 8+590436 


14015744 


14017280 


F775W,F475W 


SSTXFLS J172458. 3+591545 


14017024 


14018560 


F775W,F475W 



The IRS data were obtained in stare mode, using the 
short-low and long- low modules. The IRS short-low data 
had a background image subtracted, made by taking a 
median of all the short-low observations in the Astronom- 
ical Observation Request. The long- low data were com- 
bined by subtracting nodded pairs, then combining the 
pairs with minmax rejection. Extraction was performed 
using the SPICE tool 5 , with optimal extraction. Aper- 
ture photometry was performed on the pipeline (post- 
BCD) MIPS mosaics using standard aperture correc- 
tions 6 . 

The ACS data were taken in two filters, one longward 
and one shortward of the 4000A break in the rest frame 
of the galaxy. Observations consisted of four exposures 
in each filter in a standard dither pattern, for one orbit 
per object per filter. The ACS pipeline mosaics were also 
adequate for the analysis performed here. A full analysis 
of the data will be presented in Seymour et al. (2007, in 
preparation). 

Four of our six objects were easily classified as 
type-2 quasars in Lacy et al. (2007) based on 
well established optical emission line diagnos- 
tics. However, SSTXLFS J171106.8+590436 and 
SSTXFLS J172458.3+591525 had optical emis- 
sion lines properties which were ambiguous. For 
SSTXFLS172458.3+591525, we are able to estimate 
a ([Sii] 6719+6732) /(Ha+ [Nil]) ra tio of 0.4 from the 
IRTF/Spex spectrum, which, together with an [Om]/H/3 
ratio of 10 and assuming typical [Nll]6584/Ha ratios 
of 0.1-1.6, places this object firmly amongst the AGN 
in the diagnostic plot of Kewley et al. (2006), with a 
range in log([Sll]/Ha) « —0.4 — 0.1. In the case of 
SSTXFLS 171106.8+590436, the Ha and [Sn] lines are 
too weak to do this. However, the very high [Om]/H/3 
ratio of > 10 : 1, detection of Hell4686, and a marginal 
detection of [Nev]14.3/xm in the IRS spectrum make it 
very likely that an AGN classification for the optical 
spectrum is correct. We thus assume for the rest of this 
paper that all six objects contain powerful AGN. 

3. ANALYSIS 

We fit our SEDs and spectra using a model based on 
those of Sajina et al. (2005, 2006). The IRS spectra 
were fit directly, along with photometric points in r, ? 
and z bands from the SDSS, IRAC 3.6 and 4.5/mi data 
(Lacy et al. 2005b), MIPS 70 and 160/im data (Frayer 
et al. 2006 and this paper) and near-infrared photometry 
from IRTF. We also checked the IRS spectrophotometry 
against the IRAC 5.8 and 8.0/im and MIPS 24/im (Fadda 
et al. 2006) broad-band photometry, and found them to 
agree within w 10%. 

5 http://ssc.spitzer.caltech.edu/postbcd/spice.html 

6 http: / /ssc. spitzer.caltech.edu /mips /apercorr / 



The far-infrared (FIR) emission was fit as a modified 
black body with temperature « 45K. A warm (small 
grain) component was included with a power-law index of 
« 2, and cutoffs at high and low frequency. The FIR lu- 
minosity, I/fir, was estimated by summing the warm and 
cold components. The hot dust component was based 
on a power-law with variable spectral index. The long 
wavelength cutoff was determined by a fit of a Fermi 
function to the least starforming AGN in our sample, 
SSTXLFS J172123. 1+601214, which shows a turndown 
in the hot dust emission beginning at a rest-frame wave- 
length w 20/xm. Netzer et al. (2007) also suggest that 
the AGN-related hot dust emission starts to decrease be- 
yond ~ 20/im in type-1 quasars. The hot dust emission 
was exponentially cutoff at high frequencies with a vari- 
able cutoff frequency. The unreddened luminosity of the 
AGN, Lagn, was estimated by integrating this compo- 
nent. The hot dust component is then reddened by the 
Galactic Center extinction law of Chiar & Tielens (2006), 
likely to be a good approximation to the extreme density 
environments of AGN (Sajina et al. 2007). 

A modified version of the starburst Polycyclic Aro- 
matic Hydrocarbon (PAH) emission template of Sajina 
et al. (2007) was used to fit the PAH features. AGN are 
known to modify the PAH spectra (Smith et al. 2007). 
In particular, for our objects, we found a better match 
by boosting the 7.7/xm feature in the template by hand. 
This PAH template was fitted to the SEDs of all galaxies 
except SSTXFLS J172123.1+601214, which had a small 
amount of silicate emission added instead. Total PAH 
luminosities, Lpah, were determined by integrating over 
the fitted template. Finally, a 5Gyr old single stellar 
population from Bruzual & Chariot (2003) was fit to the 
optical points. The basic properties of the fits are listed 
in Table 2. Parameters which are not well-constrained by 
the data (for example the warm dust continuum strength 
in objects not detected at 160/im) are held fixed. Full de- 
tails of the fitting process will be given in a future paper 
(Ridgway et al. 2007, in preparation). 

4. DISCUSSION 

We first address the amount of star formation in the 
host galaxies. Star formation rates (SFRs) derived from 
FIR emission alone are suspect for quasars, as repro- 
cessed quasar emission can, in principle, contribute to 
the FIR flux. In the 5/6 cases where PAHs are detected, 
however, they show typical ratios of PAH/FIR luminosi- 
ties over a wide range of Lfir (0.02-0.06 in total PAH 
integrated over our template, corresponding to 0.0014- 
0.0042 in the 6.2/im feature, within the range for star- 
forming galaxies and ULIRGs seen by Peeters, Spoon 
& Tielens [2004] of - 0.005 - 0.1). This indicates that 
the bulk of the FIR emission is indeed likely to be from 
star formation (see also Sajina et al., in preparation). 
A similar result has been found for type-1 quasars by 
Schweitzer et al. (2006). As Figure 2 shows, however, the 
[Oil] 3727 emission (from the spectra of Lacy et al. [2007], 
aperture corrected using the SDSS photometry) is signif- 
icantly below that expected in an unreddened starburst 
in most cases, particularly for the objects with high FIR 
luminosities. This is consistent with the result of Ho 
(2005), who found low [Oil] 3727 emission from type-1 
quasar hosts compared to their expected star formation 
rates based on measurements of molecular gas. How- 
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Fig. 1. — SEDs and HST images of the six type-1 quasars, in order of increasing silicate absorption. The solid black lines are the 
IRTF/Spex and Spitzer /IRS spectra (both smoothed to reduce noise) and the black open circles photometry from SDSS, IRTF and Spitzer. 
The model fit is shown as the red dashed line, with, from left to right, the stellar component, hot dust component (including absorption, 
but excluding PAH emission) and sum of warm and cold dust components shown as blue chain-dotted lines. To the right of the SEDs we 
show expanded plots of the IRS spectra around the PAH wavelengths. The ACS images have the long wavelength image colored red and 
the short wavelength image in green. 
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Fig. 2. — The ratio of narrow line luminosity in [On] 3727 to FIR 
luminosity (circles) and in [Nell]12.8/Ltm to FIR luminosity (trian- 
gles). Host galaxies dominated by a disk component are shown as 
filled symbols, SSTXFLS J172123. 1+601214, which has an ellipti- 
cal host, is shown with open symbols. The dashed line is the ratio 
of [Oil] to FIR luminosity expected in an unextincted starburst. 
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Fig. 3. — Absorption depth in the 9.7/rni silicate feature plotted 
against axial ratio. Symbols as for Figure 2. The dashed lines 
represent the predictions of simple uniform disk models of differing 
opacity measured perpendicular to the disk. 



ever, the ratio of the [Nell]12.8//m to the FIR luminosity 
is similar for all our objects. [Nell] is of similar ioniza- 
tion to [On] and is also therefore dominated by emission 
from the starburst. This suggests that dust reddening is 
to blame for the lack of [On] 3727 emission, and is also 
consistent with the host galaxy colors and morphologies. 
As can be seen in Figure 1, the host galaxies of the ob- 
jects with high I/fir are typically disturbed, red objects, 
likely to be late stage merger remnants. This suggests 
that star formation in a significant fraction of quasar 
host galaxies is heavily obscured in the optical. 

Obscuration by dusty starforming regions also seems to 
affect the extinction to the AGN. Figure 3 shows a plot 
of axial ratio of the host galaxy (measured using the el- 



lipse task in iraf on the longer wavelength image) versus 
the depth of the 9.7/im silicate feature, S s ii, as defined 
by Levenson et al. (2006). Although our sample is small, 
the five of our six host galaxies with apparent or possible 
disk-like structures, and ld-profiles close to exponential 
disks (excluding morphological disturbances and bulge 
components), represented by the solid points, seem to 
show a trend in the sense that the more highly-extincted 
objects have more edge-on inclinations (lower axial ra- 
tios). The one galaxy which does not appear disk-like, 
SSTXFLS J172123.1+601214 (open symbol), has an un- 
derlying r 1 / 4 law profile, and is probably an example of 
an elliptical host merging with a smaller gas-rich galaxy, 
and lacks the large amount of star formation seen in the 
mergers of two gas-rich galaxies. In this case, the obscu- 
ration is probably almost entirely due to the torus. For 
the disk-like hosts, a simple uniform dusty disk model 
appears to work for most, with a spread of about a fac- 
tor of two in opacity for a line of sight perpendicular 
to the disk. One object (SSTXLFS J171106.8+590436), 
however, would have a high opacity in this model even 
if viewed face-on, corresponding to S s n ~ — 1. This may 
be because the extinction is patchy, and a high inclina- 
tion significantly increases the chance of the line of sight 
passing through a high opacity region such as a dense 
molecular cloud. Our two most obscured type-2 quasars 
(with S s [\ = —1.6 and —2) have silicate absorption as 
deep or deeper than the most obscured AGN in the sam- 
ple of Shi et al. (2006). These radio galaxies and Seyfert 
2s have S s n > —1.6, suggesting that extinction can be 
very extreme in the radio quiet quasar population. 

The implication of our result is that many type-2 
quasars have at least a component of the extinction to- 
wards the nucleus from an extended, star forming disk 
on scales of kpc, in addition to, or instead of, the tradi- 
tional obscuring torus, which interferometry suggests is 
only a few pc in size (Jaffe et al. 2004). Previous sugges- 
tions that the host galaxy can contribute significantly to 
the AGN extinction in type-2 objects have been made by 
Martinez-Sansigre et al. (2006) and Rigby et al. (2006), 
based on optical emission line properties. Our results 
suggest that a significant fraction of the quasar popula- 
tion have their AGN obscured by starforming regions in 
their host galaxies. Thus studying obscured quasars is 
essential if we are to understand the origin of the links 
between the evolution of quasars and their hosts. 
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TABLE 2 
Host galaxy properties 



Object 


z 


log(L A GN 


log(L F i R 


log(LpAH 


log(L [0 ii] 




Axial 


SFR 


SFR 






[£©]) 


[£©]) 


[Lq]) 


[Lq]) 




ratio 


(FIR) 


([Oii])t 

(Moyr- 1 ) 


















(Moyr- 1 ) 


SSTXFLS J172123. 1+601214 


0.325 


11.4 


10.6 


<9.4 


8.17 




0.63+0.02 


3.3 


3.9 


SSTXLFS J171147.4+585839 


0.800 


12.1 


11.7 


10.5 


8.44 


-0.21 


0.96+0.02 


56 


7.2 


SSTXFLS J171324. 2+585549 


0.609 


11.8 


11.3 


9.8 


8.26 


-0.34 


0.79+0.01 


16 


4.8 


SSTXFLS J171831. 5+595317 


0.700 


11.9 


12.0 


10.3 


8.45 


-0.53 


0.76+0.05 


87 


7.4 


SSTXFLS J171106. 8+590436 


0.462 


11.8 


11.5 


10.1 


7.49 


-1.56 


0.67+0.01 


33 


0.8 


SSTXFLS J172458. 3+591545 


0.494 


11.6 


11.6 


10.4 


7.91 


-1.97 


0.24+0.02 


37 


2.1 



f from Kewley et al. (2004), assuming solar metallicity. 
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